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Abstract 

The splitting of single-particle energies between spin-orbit partners in nuclei 
is examined in the framework of different self-consistent approaches, non- 
relativistic as well as relativistic. Analytical expressions of spin-orbit poten- 
tials are given for various cases. Proton spin-orbit splittings are calculated 
along some isotopic chains (O, Ca, Sn) and they are compared with exist- 
ing data. It is found that the isotopic dependence of the relativistic mean 
field predictions is similar to that of some Skyrme forces while the relativistic 
Hartree-Fock approach leads to a very different dependence due to the strong 
non-locality. 
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I. INTRODUCTION 



One often asked question about the relativistic approach to nuclear structure is: what are 
the main physical effects that a relativistic approach can handle better than a non-relativistic 
approach? Among various arguments one feature seems to emerge quite obviously, namely 
the fact that, in a description based on the Dirac equation the nucleon spin degree of freedom 
is incorporated very naturally. Therefore, nuclear properties related to the spin and par- 
ticularly spin-orbit effects in nuclei should come out from a relativistic description without 
any special adjustments whereas a non-relativistic approach would need ad hoc parameters. 
This was indeed so in the Dirac phenomenology analysis of spin observables measured in 
medium energy nucleon- nucleus scattering The interplay of two deep phenomenological 
potentials, the attractive scalar field and the repulsive vector field gives the right magnitude 
for the central and spin-orbit optical potentials. 

However, if one tries to calculate microscopically the mean field experienced by a nucleon 
bound to a nucleus and predict spin-orbit splittings, it is not so clear how successful will 
be the various relativistic models. In this work, we would like to analyze the spin-orbit 
splittings calculated in different self-consistent models: non-relativistic Hartree-Fock with 
Skyrme-type forces, relativistic mean field theory (RMFT) and relativistic Hartree-Fock 
(RHF). Of course, a mean field approach like Hartree-Fock is designed for reproducing 
global ground state properties like total binding energies or densities but it is not supposed 
to describe correctly single-particle spectra because core polarization effects are known to 
modify importantly single-particle energies PH. However, the core polarization effects 
should affect essentially in the same way the two members of a spin-orbit doublet if they 
are both below, or both above the Fermi level |3J. Therefore, it is legitimate to study the 
question of spin-orbit splittings in Hartree-Fock or mean field frameworks for nuclei where 
both spin-orbit partners are occupied. 

We carry out this study by looking at the evolution of the spin-orbit splittings of a 
nlj > — nlj < proton pair of states along an isotopic chain. The predictions of different 
models are compared with experimental values when they are available. The difficulties of 
such comparisons lie in the fact that the experimental information is scarce and it suffers 
sometimes of large uncertainties. This points to the necessity of having more data on spin- 
orbit splittings if one aims at more quantitative conclusions about the validity of different 
models concerning their spin-orbit properties. Within these limitations, our analysis will be 
done for the Ca, O and Sn isotopic chains. 

The outline of the paper is as follows. In Section II and in the Appendix we derive the 
spin-orbit potentials corresponding to the different models: Skyrme-Hartree-Fock (SHF), 
RMFT and RHF. In Section III we discuss the results for the 40 Ca- 48 Ca case as well as those 
concerning the O and Sn isotopes. Concluding remarks are made in Section IV. 

II. FORMALISM 

Let us derive the main expressions for the spin-orbit potentials in non-relativistic and 
relativistic approaches. In either case the general method is to rewrite the Hartree or Hartree- 
Fock equations for the single-particle states in a Schrodinger-like form and to identify a 
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central potential and a spin-orbit potential. These potentials are local in the cases considered 
here, but they are generally state-dependent. 



A. Non-relativistic approach 

The SHF model is very simple because the Hartree-Fock equations for the single-particle 
wave functions fa and energies take the form of a Schrodinger equation with an effective 
mass m*(r) and a local, state- independent potential. All non- locality effects are described 
by this effective mass and consequently, the SHF equations contain also first derivatives of 
the wave function. It is natural to introduce an asymptotically equivalent wave function || : 

r) , (1) 



m 

which now satisfies a Schrodinger equation with a constant mass m and a state-dependent 
local potential: 

V(e i ,r) = V (e i ,v)+V LS (r)l.s. (2) 
The spin-orbit potential is (for spherical symmetry): 

Ul* (VI rill 1 i 

V L s(q,r) = -^ -[/(r) + p q '(r)} - -[(t lXl + t 2 x 2 )J(r) + (t 2 - fi)J«(r)] , (3) 
m 1 r Ar J 

where w, ti, Xi are Skyrme force parameters, p q and J q are nucleon densities and spin 
densities (q= n or p) with p = p n + p p and J = J n + J p p. The spin densities are practically 
zero in spin-saturated nuclei but they can give some contributions in spin-unsaturated ones. 
The above expression shows that Vls is surface peaked, and that proton spin-orbit splittings 
Als should have little isotopic dependence in the SHF model as we shall see in the next 
sections. One should, however, keep in mind that another contribution to Als comes from 
the energy- dependent central potential. One can easily check that the energy dependence 
of Vo(Ei, r) is entirely contained in a term (1 — m*(r) /m)ei 0. 

Recently, some attempts have been made in order to introduce more freedom in the spin- 
orbit term of the Skyrme parametrization with the aim of improving isotope shift predictions 
in the Pb region . The authors of Ref . || modify directly the Skyrme energy functional 
so that the first term of Eq.(^) becomes wp'(r)/r + w'p q '(r)/r and the spin density terms 
are dropped. It is found that a better description of data is reached when w' is nearly — w 
(parameter set SkI4), i.e., when the proton spin-orbit potential depends almost entirely on 
the derivative of the neutron density. 



B. Relativistic approach 

Let us now calculate the spin-orbit potential in the general RHF case. The corresponding 
expressions for the RMFT case will be easily obtained by dropping all contributions corre- 
sponding to exchange (Fock) terms. Starting from the RHF equation for the 4-component 
spinors we shall obtain a Schrodinger-like equation for the upper component which allows 
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us to identify the central and spin-orbit potentials. We use the notations introduced in Ref. 
0- 

We start from Eq.(C3) of Ref. |7j for the spinor (G h Fi): 

G'i = (-- - E T>i - Pi)Gi + {m + E i + - £ ,* - 
r 

F\ = (m - Ei + + E ,i + Ri)Gi + (- + E T ,, + , (4) 

r 

where the S's are the direct (Hartree) self-energies whereas P, Q, R, S are related to the 
exchange (Fock) contributions. We shall drop the state indices i from now on. The first 
step is to eliminate the lower component F and obtain a second order differential equation 
for the upper component G: 



AG" + a x G' + a G = , (5) 



where 



A = (m + E + Zs-Zo-Q)- 1 , 
ai = A' + A(P - S) , 

a = A'(- + S T + P) + A(-4 + S' T + P') 

— (M - E + S 5 + S + i?) - A(^ + S r + P){^ + S T + S) . (6) 

Next, we look for a modified function G: 

G = \G , (7) 

such that G and G are identical asymptotically while the differential equation satisfied by 
G contains no first derivative. This condition determines A to be: 

A = CA-h-U(P-^r . (8) 

where the constant C is determined from the condition / G 2 dr = 1. In the limit of RMFT 
the exchange quantities P and S are zero and one recovers the familiar result A = A - 2. 
We can now write the differential equation satisfied by G in the form: 

The potential in state i can be identified from the coefficient of G. It is possible to separate 
out the spin-orbit part Vls of this potential by observing that the quantities A, P, Q, R, S 
depend on k = (2j + 1)(7 — j) (see, e.g., Ref. [0]) which in turn can be expressed in terms 
of (l.s). The complete expressions for Vls are given in the Appendix. 

III. EVOLUTION OF SPLITTINGS ALONG ISOTOPIC CHAINS 

We now turn to an analysis of spin-orbit splittings in some finite nuclei. We concentrate 
on the quantities: 

A(nZj< - nlj>) = e n i j< - e nW> . (10) 
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A. Ca isotopes 



In Ca isotopes, the experimental situation concerning the lc/5/2 and ld 3 /2 proton states 
is not very accurately established. In 40 Ca, the generally accepted value of ld 3 /2 single- 
particle energy is £1^3/2 = —8.3 MeV [|J, ||, [|K| However 

-15.5 MeV 



the values of £1^5/2 are much 
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E|, -16 MeV [|, -14.3 MeV pf. Note that 



more dispersed: -15.1 MeV 
the most recent values of Ref. [|1(| have been obtained as the centroid energies of the single- 
particle spectroscopic strength distributions: 



J2i EiSi 
J2i Si 



(11) 



and therefore, they take into account the strong fragmentation of the single-particle states. 
The fragment energies Ei and spectroscopic factors Si were taken from Refs. [|13|], flI5| , 
|T7p and [18| . The corresponding experimental spin-orbit splittings of 40 Ca are plotted 



on Fig. 1. In Ca, the experimental proton energies £1^3/2 are respectively -16.2 MeV 



and -15.5 MeV whereas the e ld5/2 energies are -21.5 MeV @ and -20.5 MeV [|K)f. This 
gives the spin-orbit energies in 48 Ca shown in Fig. 1. 

At this point we must mention that some theoretical papers, e.g., Ref. [JJ quote a some- 
what smaller value for the experimental A5.0. (l<i3/2 — ld5/2) in 48 Ca based on an earlier 



paper |19| but the experimental source remains uncertain. 

In Fig. [I] we display the results calculated with the following models: SHF with a 
standard SIII force and with the SkI4 parametrization of Ref. ||, RMFT with the commonly 
used NL-SH [^0| and NL3 pTJ parametrizations, RHF with (cr, lu, p, n) mesons (model e) 
of Ref. [0). The two parametrizations NL-SH and NL3 have been chosen as the most 
representative and successful in describing nuclear ground states in a wide range of nuclei. 

All calculated results for A(l<i 3 /2 — 1^5/2) show a linear A-dependence between A=40 to 
A=48. However, whereas the values of SHF and RMFT models decrease only slightly those 
of the RHF model exhibit a large decrease. 

The relatively small variations of the SIII results can be understood by examining the 
Vls potential of Eq.@. When going from 40 Ca to 48 Ca one acquires extra neutron densities 
P7/ 2 (r) and Jt/2{J') due to the filling of the I/7/2 orbital while the core neutron and proton 
densities do not change much. For the SIII parametrization, x\ = x 2 = so that J7/2O") 
gives no contribution to VlsQo, r) whereas p'j^ij) yields a contribution with a node at the 
surface, i.e., a small effect on A(ld3/2 — 1^5/2)- The same analysis holds for the SkI4 results. 
For the RMFT case the interpretation of results is also simple. The spin-orbit potential 
reduces to the first term of Eq.(|A4]), i.e., — 2A'/Ar. With the definition of A given by Eq.(P) 
it is seen that the variation of spin-orbit potential is obtained by folding the derivative of the 
scalar and vector I/7/2 densities with the (short-ranged) a and uj form factors and therefore, 
it has a node near the surface and the variation of spin-orbit splitting must be small. We 
should note, however, that the central potential that one can identify from Eq.(|9p has some 
state dependence due to the effective mass m* = m + and this effect can also contribute 
to the energy splitting. Furthermore, the kinetic energy and the central potential contribute 
also to this splitting through the j-dependence of the radial part in the wave function. 

The RHF results are more difficult to interpret. From Eq.([5|) we identify the kinetic 
energy T, central potential Vq and spin-orbit potential Vls such that: 
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(T)i + (V )i + (V LS )i = e f (l + £-) . (12) 

The expressions for Vls given in the Appendix are complicated. The non-locality of the 
original RHF mean field induces a strong state dependence in the local equivalent potentials 
and therefore, variations with j > and j < will be found in all three terms of Eq.(I2). This is 
illustrated in Table | where the values of the differences: 

AT = (T) id3/2 - (T) id5/2 , 
AV = (V ) ld3 / 2 ~ (Vq) id5/2 , 
AV LS = (V LS )ld3/2 - (Vls)u 5 /2 , (13) 

are shown. In 40 Ca the calculated spin-orbit splitting is practically equal to AVls, with a 
strong cancellation of the two other terms. In 48 Ca AVls is only one third of £1(23/2 — £1^5/2 , 
the rest coming from an increase of AVo and AT. 

It is possible to analyze further the respective role of the isovector 7r and p mesons in 
the decrease of A(lc? 3 / 2 — 1^5/2) from 40 Ca to 48 Ca. The key effect comes from the 7r meson 
and the strong non-locality related to its light mass. If one artificially increases m, up to 
600 MeV (and renormalize accordingly the f n coupling constant to keep reasonable single- 
particle levels) one tends to a more local situation and the strong decrease of A(ld 3 / 2 — lc/5/2) 
vanishes. The role of the p meson can be studied by repeating the calculations without its 
contributions to the Fock terms. It is found that the decrease of A(ld 3 / 2 — 1^5/2) still 
remains. 
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B. Oxygen isotopes 



The chain is particularly interesting because of the prospects of reaching experimen- 
tally new, unstable isotopes with a large neutron excess. We examine here the evolution 
of the lpi/2 — IP3/2 splitting of proton levels. In Fig. || we summarize the results obtained 
with SHF (SIII and SkI4 parametrizations), RMFT (NL-SH and NL3 parametrizations) and 
RHF (model e) of Ref. 0). We note that the results of NL-SH and NL3 are rather similar 
for spin-orbit splittings in isotopes. This behaviour remains also in the Sn isotopes (see 
next subsection). These two parametrizations differ mostly by their predicted compression 
modulus K and therefore, this indicates that K has little effect on spin-orbit properties in 
the RMFT. One can distinguish three intervals ending at A=22, 24 and 28 which correspond 
to the filling of lrf 5 / 2 , 2si/ 2 and subshells, respectively. The trend in the M5/2 subshell 
resembles that of the Ca isotopes, namely a relatively small variation of A(lpi/ 2 — IP3/2) for 
the SHF and RMFT models and a large decrease for the RHF model. A similar discussion as 
for the Ca isotopes can be done here. It would be interesting to measure experimentally the 
variations of spin-orbit splittings in the chain of Oxygen isotopes. For the isotopes heavier 
than A=22 RMFT shows a maximum at A=24 followed by a decrease whereas for RHF the 
behaviour is opposite, with an increase of the spin-orbit splitting after A=22. 



C. Sn isotopes 

In the chain of Sn isotopes proton spin-orbit splittings are experimentally known in 



several nuclei |22|. In Fig. y are shown the calculated and measured values of A(2pi/ 2 —2p 3 / 2 ) 
for protons, and in Fig. f| are displayed the calculated results for A(l/ 5 / 2 — 1/7/2)- From 
Fig. ^ it can be seen that none of the models is able to reproduce the very small empirical 
values of A(2pi/ 2 — 2p3/ 2 ). In the range A=112-124 the SIII results show relatively small 
variations whereas the 3 other models give large fluctuations. Of course, it would be more 
satisfactory to calculate these nuclei using a Hartree-Fock-BCS description. It is not clear, 
however, that the proton spin-orbit splittings would be significantly affected so as to bring 



them into agreement with the data. Indeed, in Ref. |23fl a HF-BCS was performed for Sn 
isotopes using the effective interaction SkI4 and a density-dependent pairing force. The 
values of A(2pi/ 2 — 2£>3/ 2 ) thus obtained vary from 1.6 MeV in 112 Sn to 2.0 MeV in 124 Sn. 
As for the I/5/2 — I/7/2 splitting, Fig. || shows a wide range of predictions. For instance, 
in A=100 RHF is as low as 2.2 MeV while the other models are in the 5-6 MeV range . It 
would be very interesting to have data on this I/5/2 — I/7/2 case i n order to shed light on 
this issue. 



IV. CONCLUSION 

We have examined in this work some predictions of spin-orbit splittings in the framework 
of non-relativistic and relativistic mean field approaches. We selected, as a representative 
case of the non-relativistic self-consistent approach, the Skyrme-Hartree-Fock model because 
it is widely used and its analytic simplicity lends itself to an easy interpretation of the 
calculated splittings. Numerical results were obtained with two versions of the effective 
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Skyrme interaction, the standard parametrization SIII and the more recently proposed SkI4 
which contains more degrees of freedom in its spin-orbit part. On the relativistic side, the 
RMFT is also a successful model for describing nuclear ground states and we have examined 
its spin-orbit predictions with the often used non-linear versions NL-SH and NL3. Here also 
the analytic form of the one-body spin-orbit potential is simple enough to allow some insight 
into the numerical results. In addition, we have included in our study the RHF model of 
Ref. [|7| in order to examine the role of the pion in the spin-orbit splittings. 

One conclusion which can be drawn is that there is no clear advantage of the RMFT over 
the non-relativistic Skyrme-Hartree-Fock approach. It is true that in RMFT, the spin-orbit 
properties come out together with the central part of the mean field whereas in the non- 
relativistic approach one has the freedom of one or more additional spin-orbit parameters. 
Nevertheless, it cannot be concluded that the RMFT spin-orbit splittings describe the data 



particularly well. In Refs. [ 24 1 , [^5| an interesting formal connection is made between RMFT 



and Skyrme-HF, but we can see here that, at a quantitative level, not only their spin-orbit 
predictions do differ, but even between SIII and SkI4 there are sizable differences in Ag.o. 

All models fail to reproduce the very small 2pi/2 — 2^3/2 proton splittings in Sn isotopes 
and it does not seem that pairing correlations can improve this prediction. As for the RHF 
approach, the light pion mass produces strongly non-local Fock fields, i.e., a strong state 
dependence of the mean fields. Consequently, the splitting between spin-orbit partners is 
due not only to (Vls^- s )ij but also to a large extent to the state dependence of the central 
potential. Thus, the evolution of the predicted splittings along an isotopic chain does not 
resemble that of less non-local models like RMFT or Skyrme. It appears that the isotopic 
dependence of A disagrees with the data in the case of RHF. 

The conclusion made earlier that RHF describes well the experimental Ag.o. both in 
40 Ca and 48 Ca was too premature in view of the uncertainties in the experimental values. 

Finally, the question of spin-orbit predictions in the framework of self-consistent theories 
is still open. Existing parametrizations of effective interactions, relativistic as well as non- 
relativistic, need further improvements. The RMF predictions seem to follow a common 
trend whereas the Skyrme-HF results may differ somewhat depending on the parametriza- 
tions, as it can be seen with SIII and SkI4 in Sn isotopes. Improving the parametrizations 
necessitates better comparisons with experiment. In this respect, it would be very helpful to 
have more data on the evolution of spin-orbit splittings along isotopic chains. Without bet- 
ter data, it is not possible at the moment to state that the spin-orbit problem is understood 
either by RMFT or by the non-relastivistic approach. 
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APPENDIX A: 



In this Appendix we give the expressions for Vls m the general case of RHF. The special 
case of RMFT is easily obtained by setting all functions P,Q,R,S to zero. The notations 
are those of Ref. |7j . 

Eq.(0) is a Schrodinger-type equation for a state of angular momentum We split 

the coefficient of G into a centrifugal term, a central potential and a spin-orbit potential: 

A" l.cti.n an —1(1 + 1) 2m rTr/ , . I-, ,. . 

T " 2 ( T } + X = r 2 " "^rtW + Vbs(r)a.s)y - e(l + e/2m)]. (Al) 

The functions P,Q,R,S can be expressed as second order polynomials in k , namely, P = 
P + Pxk + P 2 k 2 and similarly for Q, R, S. Here k = (2j + l)(l - j) = |Q(2j + 1) where 
D, — +1 if j — I — | and Q — — 1 if j — I + |. It is easy to relate powers of k to (l-s)y. For 
instance, we have: 

2(l.s) y = -h 2 (l + k ). (A2) 

Thus, all contributions to Vis(r) can be identified by expressing the K-dependence into a 
(l.s)jj-dependence. One then obtains: 

Vls{t) = vg + + , (A3) 

where 

vi3 = I{-£(i + *) + ^-j* 1 + * + (; + ^)(^ + ^) + (; + 5 1 )(e- + p ) 



+^P 2 + P' 2 - ^ - (5# + Po)S 2 - (~ + Pi) (~ + Si) - P 2 (£? + So) 

+ [l(l + !) + !} [(- + Pi)S 2 + (- + S^P,] - [2/(/ + 1) + 1]P 2 S 2 } , (A4) 



= ±{- ^(P - Sx) - (P - S )(P a - Si) 

+^(P 2 - S 2 ) + i(Pi - Si) 2 + (P - S )(P 2 - S 2 ) 

-[/(/ + 1) + 1](P - Sx)(P 2 - S 2 ) + [2Z(/ + 1) + l]i(P 2 - S 2 ) 2 } , (A5) 
^{ - ^(P ~ Si) + ^(P'i - S\) - ^(P - S )(P 1 - Si) 

- S 2 ) - ~(P' 2 - S' 2 ) + ~(P X - Si) 2 + ~(P - S )(P 2 - S 2 ) 

-[/(/ + 1) + l]i(Px - Si)(P 2 - S 2 ) + [2Z(Z + 1) + 1]-{P 2 - S 2 ) 2 } . (A6) 



y(UI) 
V LS 
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TABLES 



TABLE I. The quantities AT, AVb, AVls and £^3/2 — £\d5/2 calculated with RHF for protons 
in 40 Ca and 48 Ca (in MeV). 



Nucleus 


AT 


AL 


AV LS 


£l(23/2 — £ld5/2 


40 Ca 


-2.94 


3.31 


7.54 


8.01 


48 Ca 


-1.73 


4.25 


1.38 


4.06 
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FIGURES 



FIG. 1. The A(lt2 3 / 2 — W5/2) proton splittings in Ca isotopes calculated with SIII (dot-dashed 
line), SkI4 (long-dashed line), RMFT (dotted line) and RHF (full line) models. Experimental values 
in A=40 and A=48 are shown: A [10], B [11], C [8] and D [9]. 

FIG. 2. The A(lp 1 / 2 — IP3/2) proton splittings in O isotopes calculated with SIII (dot-dashed 
line), SkI4 (long-dashed line), RMFT (dotted line) and RHF (full line) models. The experimental 
value for A=16 is indicated. 

FIG. 3. The A(2p 1 / 2 — 2^3/2) proton splittings in Sn isotopes calculated with SIII (dot-dashed 
line), SkI4 (long-dashed line), RMFT (dotted line) and RHF (full line) models. Experimental 
values for A >112 are shown. 

FIG. 4. The A(l/ 5 / 2 — I/7/2) proton splittings in Sn isotopes calculated with SIII (dot-dashed 
line), SkI4 (long-dashed line), RMFT (dotted line) and RHF (full line) models. 



12 



l — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 




j i i I i i i i I i i i i I i i i i I i i i L 



Z 9'L I S'O 

( A 8fl) ( z/£ d Z - z/l d Z )v 



